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Abstract. We present spatially resolved Chandra HETGS observations of the Seyfert 2 galaxy NGC 1068. X- 
ray imaging and high resolution spectroscopy are used to test the Seyfert unification theory. Fe Ka emission 
is concentrated in the nuclear region, as are neutral and ionized continuum reflection. This is consistent with 
reprocessing of emission from a luminous, hidden X-ray source by the obscuring molecular torus and X-ray 
narrow-line region (NLR). We detect extended hard X-ray emission surrounding the X-ray peak in the nuclear 
region, which may come from the outer portion of the torus. Detailed modeling of the spectrum of the X-ray NLR 
confirms that it is excited by photoionization and photoexcitation from the hidden X-ray source. K-shell emission 
lines from a large range of ionization states of H-like and He-like N, O, Ne, Mg, Al, Si, S, and Fe xvii-xxiv L-shell 
emission lines are modeled. The emission measure distribution indicates roughly equal masses at all observed 
ionization levels in the range log^ = 1—3. We separately analyze the spectrum of an off-nuclear cloud. We find 
that it has a lower column density than the nuclear region, and is also photoionized. The nuclear X-ray NLR 
column density, optical depth, outflow velocity, and electron temperature are all consistent with values predicted 
by optical spectropolarimetry for the region which provides a scattered view of the hidden Seyfert 1 nucleus. 
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1. Introduction 

NGC 1068 (z=0.003793) contains one of the brightest and 


comes from fluorescence in the Compton-thick torus, with 
a possible contribution from the ionized X-ray narrow- 
line region (NLR) (Krolik fc Kallman 1987 ; Marshall et 


served as the e xemplar of AGN uiiiricaLiu ii fur Seyfert 1 
and 2 galaxies ( Antonucci fc Miller 1985 ). Broad optical 
lines from a hidden Seyfert 1 nucleus are scattered by elec¬ 
trons in the nuclear region. It is hypothesized that Seyfert 
1 and 2 galaxies are primarily distinguished by whether 
we view the nucleus directly or indirectly by scattered 
and reprocessed light. This model can be tested in many 
ways using high resolution spatial and spectral observa¬ 
tions with the Chandra X-ray Observatory. 

One of the predictions of the Unified Seyfert Model is 
the existence of an optically thick torus which collimates 
the ionizing radiation from the nucleus and blocks our di¬ 
rect view of the nuclei of Seyfert 2 galaxies. Observations 
of the hard (E > 10 keV) X-ray continuum of NGC 1068 
with Beppo-SAX show that it is dominated by Compton 
r eflection, and sug gest that the torus is Compton-thick 


closest|3eyfeil 2 active galacHc imdC (AGN). lUias also |[I^ ^att et ah 1119961). Our high-resolution Chandra 


alent width Fe Ka emission line ( 

Koyama et al 

1989; 

Marshall et al. 

1993; 

Iwasawa et al. 

1997) most likely 
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observations allow us to constrain the location, size, and 
kinematics of the Fe Ka emitting region and investigate 
the properties of the torus. 

It has been suggested that warm ionized gas ablated 
from the torus provides the optical scattering medium 


(Krolik & Begelman 1986; Krolik & Kriss 1995). One 


important implication of electron scattering by warm gas 
in the NLR is that there must be a phase of the NLR 
with temperature T ~ 10®K. This plasma should also re¬ 
veal itself in the X-ray band via scattered continuum, re¬ 
combination emission, and resonance scattering. We make 
high-resolution measurements of these spectral compo¬ 
nents with Chandra HETGS, and derive strong constraints 
on the properties of the electron scattering region. 

Another important goal is to understand the physi¬ 
cal state and dynamics of AGN outflows. The bi-conical 
narrow-line region (NLR) seen in many Seyfert galax¬ 
ies is one manifestation of these outflows. Optical spec¬ 
troscopy of the NGC 1068 NLR with the Hubble Space 
Telescope {HST) has revealed the kinematics of the inner 
NLR, which appears to be an outflowing wind (prensh^ 
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& Kraemer 2000). Interaction of this wind with the radio 
jet and h ost galaxy ISM may also be important ( Cecil el 
al. 


2001 ). 


Recent Chandra imaging spectroscopy of NGC 1068 
shows a detailed correspondence between extended X-ray 


emission and the optical narrow-line region (Young et al. 


2001) Similar extended X-ray nebular emission is ob¬ 


served in NGC 4151, Mrk 3, and Circinus (pgle et al. 


200C; Sako et al. 2000| ; Sambruna et al. 1 |2001| ). XMM 
RGS observations of the soft X-ray spectrum of NGC 1068 
show the importance of photoexcitation and give the first 
measurements of the temperature and column density of 
X-ray emitting gas in the NLR (Kinkhabwala et al. 2002] ). 
We apply a similar analysis to Chandra HETGS spectra 
of NGC 1068, which have better resolution and S/N in the 
medium-hard (0.8-7 keV) band. 

The high spatial resolution of Chandra HETGS en¬ 
ables us to measure the spectrum in two locations in the 
NLR-the central 1.5"of the nuclear region, and a cloud 
3" NE of the nucleus (subsequently, “NE cloud”). We use 
these spectra to measure and compare the column den¬ 
sities in the two regions. Multi-band X-ray imaging also 
reveals how the X-ray spectrum changes across the NLR. 
Near-simultaneous observations with Chandra LETGS are 


presented in a parallel paper by Brinkman et al. (2002). 


2. Observations 

We observed NGG 1068 on 2000 December 4-5 (MJD 
51882-3) with Chandra HETGS. The total exposure time, 
after correcting for dead time is 47.0 ks. The disper¬ 
sion axes (Roll= 308.9°, HEG PA= 123.7°, MEG PA= 
133.6°) are nearly perpendicular to the axis of the NLR 
(PA~ 40°). This maximizes the spectral resolution and 
enables spectroscopy at multiple locations along the ex¬ 
tended NLR. The width of the nuclear emission region in 
the dispersion direction for the zeroth order image ranges 
from FWHM= 0.81 - 0.66" for A = 0.6 - 3 keV. This cor¬ 
responds to a smearing of FWHM = 0.015 — 0.018 A over 
the 6-22 A range in addition to the instrumental profile 
(FWHM= 0.01, 0.02 A for HEG, MEG). 

Events were filtered by grade (eliminating grades 1, 
5 and 7) and streak events were removed from chip S4. 
First-order HEG and MEG spectra of the nuclear region 
(Fig. 0) and off-nuclear region (Fig. ||) were extracted us¬ 
ing GIAO 2.1. Extraction windows in the cross-dispersion 
direction are 3.2" and 3.0", respectively. The nuclear ex¬ 
traction window is centered on the peak X-ray emission. 
The off-nuclear extraction window is centered on a cloud 
3.1" NE of the nucleus (NE cloud). The two spectral ex¬ 
traction regions are adjacent and do not overlap. Spectra 
were divided by grating effective area curves generated 
from the aspect histogram and instrument effective area. 

The ACIS zeroth order image provides useful informa¬ 
tion about the spatial and spectral distribution of X-ray 
emission. The central 3x3 pixels have a count rate of 
0.31 counts per 3.2 sec frame, which indicates moderately 
heavy pileup. This precludes accurate spectral analysis of 


the X-ray peak, but should not greatly affect the surround¬ 
ing regions. We have split the zeroth order image into 6 
energy bands: 0.4-0.6, 0.6-0.8, 0.8-1.3, 1.3-3, 3-6, and 6-8 
keV (Fig. ^). These bands correspond roughly to the O 
VII, O VIII, Ne & Fe L, Mg Xi-S xv, hard continuum, and 
Fe K line series, respectively. The individual images are 
30" = 2.2h^^ kpc on a side. We adaptively smoothed the 
broad-band images using the GIAO program GSMOOTH 
and combined them into two 3-color images (Figs. 0, ||), 
which are 20" = 1.5h^^ kpc on a side. The high-energy 
image (Fig. 0) shows the distribution of Fe Ka, scattered 
continuum, and high-ionization line emission. The low- 
energy image (Fig. 0) shows the distributions of lower 
ionization states and Fe L emission. 


3. Nuclear Spectrum 

The nuclear spectrum of NGC 1068 (Fig. 0) consists of 
emission from H-like and He-like ions of C, N, O, Ne, Mg, 
Al, Si, S, and Fe, and Fe xvii-xxiv. Lines were fitted with 
Gaussians to determine their rest wavelengths and fluxes 
(Table 0 , ||). Positive IDs were made for lines within 0.04 
A of the predicted wavelengths for abundant ions. 

We fit the brightest unblended emission lines from the 
MEG spectrum of the nuclear region with Gaussian pro¬ 
files to determine their redshifts and Doppler b parame¬ 
ters (Table ||) . The profiles are convolved with a Gaussian 
spatial profile, with FWHM= 0.7". We find that all of the 
lines are blue-shifted, with < v >= —410 ± 80 km s“^. 
This is consistent with the HST [O ill] spectrum of the 
nuclear region ( prenshaw fc Kraemer ^000 ), and shows 
that the X-ray emitting plasma is flowing out along with 
the optical NLR clouds. All of the H-like and He-like X-ray 
lines from the nuclear region are resolved. We find a range 
of line widths, from b = 690 km s“^ for O vii to 6 = 1600 
km s“^ for Mg xii (Table 0 ). These widths are signif¬ 
icantly larger than the average optical line widths, but 
consistent with the most extr eme [O ill] veloci ty widths 
found in optical NLR clouds ( pecil et al. 2001 ). 

Photoionization is indicated by the strength of the for¬ 
bidden lines and prominent narrow radiative recombina¬ 
tion continua (RRCs). The ratios r/f of resonance to for¬ 
bidden lines from the He-like ions (e.g. Mg xi 2p-ls/f) are 
intermediate between those ex pected for pure recombina - 
tion and collisional excitation ( Porquet fc Dubau ^000| ). 
The nuclear spectrum of NGC 1068 is very similar in 


this r espect to the spectrum of NGC 4151 ( Ogle et al 
|2000| ). However, a large value of r// doesn’t by itself in¬ 
dicate a collisionally ionized component to the spectrum. 


As demonstrated with XMM RGS (Kinkhabwala et al. 
[2002 ), photoexcitation enhances the high-order resonance 
transitions as well as the n = 2 — 1 (r) line. The strong 
high order transitions in the Chandra HETGS spectrum of 
NGC 1068 (Fig. 0 , Table 0 , 0 ) confirm that photoexcita¬ 
tion, not collisional excitation is responsible for the strong 
n = 2 - 1 (r) lines of Mg XI and Si xiii. 

Fluorescence is also an important contributor to the 
nuclear spectrum. We see strong, narrow Fe i Ka and de- 
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Energy (keV) 

8 7 6 5 4 3 



2.2 2 1.8 1.6 1.4 1.2 1 



0.8 0.75 0.7 0.65 0.6 0.55 0.5 



Wavelength (A) 

Fig. 1. Chandra HETGS spectra of the central 1.5" radius nuclear region. The top panel is from the HEG (0.005 A 
bins); bottom two are from MEG (0.01 A bins). The strong forbidden emission lines (f) and narrow recombination 
continuum features (RRG) indicate recombination following photoionization by the hidden nucleus. Fe I, Si I, and S 
I Ka emission lines indicate fluorescence from a low-ionization region, which may be identified with the obscuring 
molecular torus. 


tect Fe I K/3, S i Ka, and Si i Ka lines. Fe i Ka is ob¬ 
served at a rest wavelength of A = 1.934 ± 0.02 A , and 
the core is unresolved with a width of cr = ^ 10“^ 

A . If the line is purely from Fe i, it is significantly blue- 
shifted, with V = —500 ± 300 km s“^ and its width corre¬ 


sponds to a turbulent velocity of 6 < 1100 km s“^ (90% 
confidence). Alternatively, the line may be a blend of Fe 
I-XVIII (Decaux et al. 1995), which could cause an appar¬ 
ent blue-shift. The S i Ka line is also blue-shifted, with 
V = —390 ± 170 km s“^. Fe i Ka has a weak red Compton 
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Energy (keV) 

8 7 6 5 4 3 



2.2 2 1.8 1.6 1.4 1.2 1 



0.8 0.75 0.7 0.65 0.6 0.55 0.5 



Wavelength (A) 

Fig. 2. Chandra HETGS spectra of the NE cloud, centered 3.1" NE of the nucleus. Binning is the same as in Fig. 1. 
The forbidden lines (f) and RRCs are weaker than in the nuclear region (Fig. 1), due to lower column density. The 
H-like and He-like resonance lines are relatively stronger from photoexcitation. The Fe L emission lines are stronger, 
indicating a greater Fe abundance. Fe Ka emission is weaker because there is relatively little neutral matter in this 
region. 


shoulder, as observed by Iwasawa et al. (1997). The ra¬ 


tio of Fe K/3 to Fe Ka is 0.17 ± 0.06, consistent with the 
predicted ratio (~ 0.13) of fluorescence yields (Kaastra & 
Mewe 


1993). 


The hard nuclear continuum emission (2.5-6.2 keV) is 
roughly fit by a power law with T = 1.2 ± 0.3 and F 2-10 = 
2.8 X 10-12 erg cm-^ s"! (L 2-10 = 7.7 x lO^o erg s"!), 
using a (fixed) Galactic column of A^h = 3.5 x lO^^ cm-^. 
This is harder than the direct spectrum (T = 1.8 ± 0.2) 
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Table 1. NGC 1068 emission line fluxes (10 ® ph s ^ 
cm“^) from the nuclear region. Predicted wavelengths 


(|1985|),[Drake ( 

1988), Brown et al. 

(1998), and Brown et 

al. (|2C|02). Observed wavelength is 

divided by 1 -|- Zg^i. 

Ap 

-^obs 

± 

Flux 

± 

ID 

1.757 

1.755 

0.005 

1.2 

0.5 

Fe I K/1 

1.780 

1.797 

0.007 

0.9 

0.6 

Fe XXVI Lya 

1.850 

1.843 

0.005 

1.4 

0.6 

Fe XXV 2p-ls (r)** 

1.858 

1.854 

0.012 

1.2 

0.7 

Fe XXV i ® 

1.868 

1.881 

0.005 

2.1 

0.6 

Fe XXV f 

1.9373 

1.9357 

0.0010 

5.1 

0.6 

Fe I Ka 

5.039 

5.045 

0.009 

0.6 

0.2 

S XV 2p-ls (r) 

5.086, 

5.101 

5.099 

0.005 

0.5 

0.2 

Si XIII RRC -t S 

XV f 

5.374, 

5.405 

5.367 

0.003 

0.4 

0.2 

S I Ka -1- Si XIII 
4p-ls 

5.681 

5.676 

0.004 

0.3 

0.1 

Si XIII 3p-ls 

6.182 

6.181 

0.002 

1.0 

0.1 

Si XIV Lya 

6.648 

6.636 

0.002 

1.2 

0.2 

Si XIII 2p-ls (r) 

6.69 

6.676 

0.004 

0.6 

0.1 

Si XIII i ==■ 

6.740 

6.728 

0.002 

1.1 

0.1 

Si XIII f “ 

7.106 

7.096 

0.006 

0.6 

0.1 

Mg XII LyP 

7.130 

7.128 

0.004 

0.6 

0.1 

Si I Ka 

7.310 

7.318 

0.005 

0.04 

0.04 

Mg XI 5p-ls 

7.473 

7.480 

0.015 

0.2 

0.1 

Mg XI 4p-ls 

7.757 

7.746 

0.003 

0.20 

0.07 

Al XII 2p-ls 

7.850 

7.839 

0.005 

0.5 

0.1 

Mg XI 3p-ls 

8.421 

8.397 

0.005 

2.8 

0.3 

Mg XII Lya 

9.169 

9.159 

0.004 

1.7 

0.2 

Mg XI 2p-ls (r) 

9.23 

9.210 

0.008 

0.7 

0.2 

Mg XI i** 

9.314 

9.297 

0.004 

1.2 

0.2 

Mg XI P 

10.007 

10.006 

0.007 

0.5 

0.1 

Fe XX 


De-blended 


from a typical Seyfert 1 galaxy (Risaliti 

20021), and in- 

dicates neutral Compton reflection, consistent with the 

strong Fe i Ka fluorescence line (Matt et al. 

19961). The 


reflection spectrum is analyzed in more detail below. 


4. Off-Nuclear Spectrum 

The spectrum of the NE cloud (Fig. ||) is quite different 
from the nuclear spectrum (Fig. One striking differ¬ 
ence is the weakness of the Fe i Ka line. This indicates 
that most of the fluorescence occurs in the nuclear region 
and not in the extended NLR. 

Recombination emission is also considerably weaker in 
the NE cloud than in the nuclear region. Note the weak¬ 
ness of the O VII f and Ne ix f forbidden lines and the O 
VII RRC in Fig. || This indicates a relatively low column 
density of plasma in the NE cloud, which we model in de¬ 
tail below. Emission from Fe xvii-xxiv L lines dominates 
the 0.7-1.3 keV band. Strong Fe L lines, (e.g. 3d-2p lines 
of Fe xvii) could be mistaken as an indicator of collisional 
excitation. However, as we show below, these lines arise 


Table 2. Continuation of Table 1. 


Ap 

Aobs 

± 

Flux 

± 

ID 

10.239 

10.228 

0.005 

1.3 

0.2 

Ne X Ly/3 " 

11.000, 

10.984 

0.004 

1.3 

0.2 

Ne IX 4p-ls 

10.981 





-1- Fe XXIII 3p-2s 

11.547, 

11.520 

0.004 

1.2 

0.2 

Ne IX 3p-ls 

11.490 





-I- Fe XXII 3p-2p 

12.134 

12.114 

0.004 

3.5 

0.4 

Ne X Lya 

12.284 

12.257 

0.004 

0.9 

0.2 

Fe XXI 3d-2p 

13.447 

13.430 

0.006 

2.2 

0.6 

Ne IX 2p-ls (r) 

13.518, 

13.520 

0.006 

2.8 

0.5 

Fe XIX 3d-2p -1- Ne 

13.55 





IX i 

13.698 

13.680 

0.003 

6.0 

0.6 

Ne IX f 

13.795 

13.777 

0.005 

0.8 

0.3 

Fe XIX 3d-2p 

14.034 

14.06 

0.05 

0.7 

0.6 

Fe XIX 3d-2p 

14.208 

14.178 

0.008 

4.4 

0.7 

Fe XVI 11 3d-2p 

14.373 

14.36 

0.02 

1.4 

0.6 

Fe XVI 11 3d-2p 

14.534 

14.51 

0.02 

1.9 

0.7 

Fe xviii 3d-2p 

14.961 

14.963 

0.005 

1.6 

0.5 

Fe XIX 3s-2p 

15.014 

15.024 

0.011 

1.4 

0.5 

Fe XVII 3d-2p 

15.079 

15.053 

0.012 

0.5 

0.5 

Fe XIX 3s-2p 

15.176 

15.160 

0.011 

1.2 

0.5 

O VIII Lyy 

15.261 

15.232 

0.002 

1.3 

0.3 

Fe XVII 3d-2p 

15.453 

15.447 

0.009 

0.3 

0.2 

Fe XVII 3d-2p 

16.006 

15.975 

0.005 

1.0 

0.4 

O VIII Ly/3 

17.051 

17.035 

0.017 

3.4 

1.0 

Fe XVII 3s-2p 

17.200 

17.161 

0.012 

0.6 

0.4 

O VII 6p-ls 

17.396 

17.376 

0.006 

0.3 

0.4 

O VII 5p-ls 

17.768 

17.758 

0.004 

2.3 

0.7 

O VII 4p-ls 

18.629 

18.60 

0.02 

1.4 

1.2 

O VII 3p-ls -1- N 

VII RRC*" 

18.969 

18.947 

0.014 

7.7 

1.9 

O VIII Lya 

20.910 

20.883 

0.008 

2.6 

1.2 

N VII Ly/3 

21.602 

21.589 

0.007 

8.5 

2.1 

O VII 2p-ls (r) 

21.804 

21.747 

0.005 

1.8 

1.3 

O VII i 

22.097 

22.072 

0.002 

20.4 

3.5 

O VII f 

24.781 

24.741 

0.008 

8.0 

2.1 

N VII Lya 


** De-blended 


from recombination and photoexcitation, as do the lines 
of the H-like and He-like ions. 

We fit unblended lines from the NE cloud with 
Gaussians, convolved with a Gaussian spatial profile with 
FWHM= 1.5" (estimated from the 0.7-1.3 keV image). 
There is no measurable redshift, and the mean velocity is 
< V >= 0 ± 330 km s“^ with respect to the host galaxy 
rest frame. This is consistent with the low (red-shifted) 
velocities meas ured for optical [Om] e mitting clouds in 
this region by Crenshaw & Kraemer ( 2000]) . They in¬ 
terpret this as evidence for deceleration by the ambient 
medium in the host galaxy, starting at a distance of 1.7" 
from the nucleus. Alternatively, this cloud may mark the 
shock front where the radio lobe impacts the ISM of the 
host galaxy (papetti et al. 1997; Pecontal et al. 1997; 
pecil et al. pOOl ). Three emission lines are resolved (O 
VIII Lya, O vii f, and N vii Lya), with < b >= 810 ± 30 
km s“^, similar to values in the nuclear region. A fourth 
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Table 3. NGC 1068 Nuclear Region. Ion velocity v and 
Doppler parameter b (km s“^) are from a Gaussian fit. 
Electron temperature Te(10'^ K) is from an RRC fit. 
Predicted temperature T^ p is from an XSTAR photoion¬ 
ization model. 


Ion 


Line 


C VI 





2.8 

N VI 





3.2 

N VII 

Lya 

-350i“° 

10401111 

Q+OO 

5.1 

O VII 

f 

-530111“ 

600111“ 

81“ 

3.2 

0 VIII 

Lya 


12301111 


8.5 

Ne IX 

f 

-370111“ 

7501111 

9l“ 

8.5 

Ne X 

Lya 

-46011““ 



20 

Mg XI 

f 

-3001111 

13801111 

y+20 

' -4 

13 

Mg XII 

Lya 

-4401111 

15701111 


33 

Si XIII 

f 

-5401111 

6901111 


33 

Si XIV 

Lya 

-3301111 

740l«l 


50 


(Wilson & Ulvestad 1987; Capetti et al. 1997; Young el 
al. I 20011 ). The bright X-ray cloud 3" NE of the nucleus 
(NE cloud) lies in the back-flow region of the radio lobe, 
but does not correspond to any individual cloud in the Ha 
or [O III] images. The NE cloud dominates the off-nuclear 
HETGS spectrum (Fig. ||). 

The spiral structure of the Ha emission appears to be 
associated with the central few arc seconds of the starburst 


disk in the host galaxy ( [Capetti et al. 1997). The NE cone 
of the X-ray NLR is seen projected in front of the disk, 
while the SW cone is hidden behind it. The SE spiral arm 
in the 0.8-1.3 keV X-ray image (Fig. ||) matches the arm 
seen in Ha, consistent with energetic starburst activity. 
The starburst disk also extends out to much larger radius, 
which is much better seen in the ACIS image of Young et 

~ar|( p00l| ). 


5.2. Hard X-ray Image 

The image in the 6-8 keV band (Fig. ||) is dominated by 
Fe I Ka emission, though there is also significant contribu¬ 
tion from Fe xxv (see Fig. [i|). This emission is spatially 


line, Mg Xii Lya, has a rather large Doppler parameter of 
b = 2900 ±500 km s“^, but may be blended with Fe xxiii. 


5. Broad Band Images 

5.1. Comparison to Ha 

To put our X-ray images in the context of the vast lore 
of optical emission line studies of NGC 1068, we have 
combined them with an archival HST WFPC2 Ha image 
(Fig. 1^). The HST image was taken through the F658N 
filter for 900 s, and the e stimated continuum (F791W fil¬ 
ter) was subtracted out ([Capetti et al. 1997; Kishimoto 


resolved, as first shown by [Young et al. | ( 2001 ). There is 
a weak (~ 10") linear extension in the 6-8 keV band (Fig. 
^) corresponding to the extended ionization cones, seen 
on both sides of the nucleus. The peak of Fe K emission is 
close to the hidden nucleus (Fig. p, and may correspond 
to the illuminated inner wall of the torus. 

The X-ray peak is surrounded by extended (2.7" x 
2.2") emission in the 3-8 keV band (Fig. H). We attribute 
this to the outer edge of an irregular molecular torus seen 
in scattered hard X-rays. This feature corresponds to a 
region of high optical extinction in HST images and lies 
at the sa me radius (~ 1") as the m olecular ring seen in CO 
emission (jSchinnerer et al~ 2000| ). The location and size 
of the hard X-ray ring is also consistent with shadowing 
of the far (SW) scattering cone observed in the near-IR 


and attributed to the torus (Young et al. 1996). 


1999a ). Red represents Ha ± N ii, green soft X-rays (0.4- 
0.7 keV), and blue the 0.7-1.3 keV band. The p eak of 


the X-ray im age was aligned to match cloud B ( Evans 
et al. I 1991 ) in the Ha image. This is the closest cloud 
to the hidden nucleus, as determined by UV polarime- 


Emission from the ionization cones is also seen in the 
1.3-3 keV and 3-6 keV bands (Fig. |3[ green and red in Fig. 

The 1.3-3 keV image is dominated by emission from 
highly ionized Mg, Si, and S, while the 3-6 keV emission is 
mostly scattered continuum. The 3-6 keV continuum pho¬ 
tons (green) penetrate the ISM of the host galaxy disk, 
revealing the far-side ionization cone which is mostly in¬ 
visible at lower energies. This is consistent with the Wh 


try ( Kishimoto 1999b). N ote that this is slightly different 
from Young et al. ( 2001| ) who align the Chandra ACIS 
images with the radio core. However, this makes little dif¬ 
ference at Chandra resolution (0.5"). 

Chandra ACIS images ( Young et al. ^001 ) show a 
strong correspondence between the optical [O HI] and X- 
ray emission in the NLR. We find a weaker correspondence 
between the Ha and X-ray emission. In fact, it appears 
that the highly ionized X-ray NLR is surrounded by a 
cocoon of Ha filaments. The tip of this cocoon, 5" NE 
of the nucleus, corresponds to the tip of the radio lobe 


map produced by Young et al. 

to 

o 

o 

1— > 

, and near-IR imag- 

ing which first revealed the SW 

cone 

[Young et al. 

1996 


Packham et al. 1997). The far cone is mostly obscured at 
optical-UV and soft X-ray energies, though the cloud 1.5" 
SW of the nucleus (Fig. ||) peeks through. 

5.3. Soft X-ray Image 

The X-ray color varies significantly across the low-energy 
image (Fig. H), owing to variations in emission line 
strengths. These variations can be attributed to ioniza¬ 
tion or abundance variations. The region extending 1.5" 
NE of the nucleus is pink, owing to relatively weak 0.8-1.3 
keV emission, and suggesting a relatively low ionization 
parameter. This corresponds closely to the strongest [O 
III] emission (arrowhead sha ped) from the nuclear region 
in 7757'images ( Evans et al. [1991 ; Macchetto et al. 1994 ). 

The NE cloud changes in color from blue to green to 
red along the radius vector from the nucleus (Fig. H), 
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Fig. 3. Six broad-band images filtered from the HETGS 
zeroth order image and adaptively smoothed. First row: 6- 
8, 3-6, 1.3-3 keV. Second row: 0.8-1.3, 0.6-0.8, 0.4-0.6 keV. 
Each image is 30" = 2.2h7g^kpc on a side, 0.49"/pixel. 
Note the change in structure with energy. Fe Ka emis¬ 
sion is concentrated in the nuclear region in the 6-8 keV 
image. Subtle differences between the various soft X-ray 
bands indicate different spatial distributions for various 
ionization states. Starburst emission from the host galaxy 
disk is most apparent east of the nucleus in the 0.8-1.3 
keV image. 


perhaps indicating a change in the ionization level. Red 
regions in Fig. ^ correspond to regions of strong Ha emis¬ 
sion (Fig. which is indeed consistent with lower ioniza¬ 
tion. The NLR clouds 1.5" SW, and 1" W of the nucleus 
(Fig. I) show a similar effect. The lowest ionization (red) 
regions in the NE, SW, and W clouds may indicate the in¬ 
tersection of the ionization cone with high density clouds 


in the disk of the host galaxy (See Fig. 6 of Cecil et al. 

(S))- 

The blue gap between the nuclear region and the NE 
cloud is filled with strong 1 keV emission. There is also 
strong Ha emission in this region (Fig. ^). There are a 
few possible interpretations. 1) This could be a region of 
low density and high ionization which emits strongly in 
Ne X and Fe L lines. If so, the Ha filaments which have 
lower ionization may just be seen in projection. 2) The 
emission could come from collisionally ionized plasma in 
the starburst disk. The Ha emission could then be from 
star forming regions. 3) This region could be shock-heated 
by the NLR outflow or radio jet. The radio jet appears to 
be deflected in this region, which m ay lend support to this 
interpretation ( Young et al. 2001 ). 


6. Discussion 

6.1. Photoionization Model 


As seen in other Seyfert 2s ( ^ako et ^ 200C| ; pambruna el 
al. 2001), photoionization, recombination, and photoexci- 



Fig. 4. Three-color high energy X-ray image {E =1.3-3, 
3-6, 6-8 keV)= (r,g,b). The image is 20" = kpc on 

a side. The brightest point-like source of high-energy emis¬ 
sion may be the inner wall of the molecular torus, reflect¬ 
ing X-rays from the hidden nucleus. The hard emission 
surrounding the hard X-ray peak may delineate reflection 
from the outer reaches of the obscuring torus. Extended 
emission to the NE and SW corresponds to the ionization 
cone. The SW cone is seen through and partly obscured 
by the disk of the host galaxy. 


tation appear to dominate the spectrum of NGC 1068. The 
main indicators are: (1) narrow recombination continua of 
N VII, O VII, and Ne ix, (2) strong forbidden lines of O vii 
and Ne ix (3) strong np-ls (n>2) lines from O vii, O viii, 
Ne IX, Mg XI, and Si xiii, (4) weak Fe L emission (rela¬ 
tive to collisionally ionized plasmas). We therefore start 
with a pure photoionization model. In a subsequent sec¬ 
tion we assess the possibility of additional emission from 
a collisionally ionized component. 

We model the nuclear and off-nuclear spectra of NGC 
1068 using the method of Kinkhabwala et al. (^002[). The 


geometry is a collection of clouds photoionized by an ex¬ 
ternal X-ray source. Each cloud is characterized by its ra¬ 
dial velocity, turbulent velocity, temperature, column den¬ 
sity, and source covering fraction. The source is character¬ 
ized by its 2-10 keV X-ray luminosity Lx and power-law 
spectral index, which we fix at T = 1.8, a typical value for 
Seyferts. The velocity and turbulent velocity determine 
the line redshifts and widths. The electron temperature 
Te controls the width of the radiative recombination con¬ 
tinuum. The column density Ni and Doppler b parameter 
determine the relative line contributions from recombina¬ 
tion and photoexcitation. Note that the photoexcitation 
rate saturates at a larger column density for larger val- 
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Fig. 5. Three-color low energy X-ray image {E =0.4-0.6, 
0.6-0.8, 0.8-1.3 keV)=(r,g,b), showing the distribution of 
line emission, which peaks near the nucleus. The color of 
the X-ray NLR ranges from blue to red, corresponding to 
high through lower ionization states. Same scale as Fig. 


ues of b. Line intensity is proportional to column density, 
covering fraction /(., and X-ray luminosity. 

We fit the spectra to derive the flux-weighted mean 
cloud parameters for each ion. Spectral fitting is done 
with IMP, a code developed at UCSB which employs a 
minimization algorithm. For H-like and He-like ions, 
(n = 1 — 7, 1) level populations are computed by solving 
rate equations including terms from photoionization, re¬ 
combination, and photoexcitation ( Behar et al. 2001 ). 
Recombination rates for n > 7 are extrapolated using 
the hydrogenic approximation. The emission line spec¬ 
trum of each ion is computed from the level populations 
and transition rates. The radiative recombination con- 


tinua (RRCs) are computed using the method of Liedahl 
( 1999[ and the photoionization cross sections of Verner 


et al. Tl|l9^) . Atomic data, including energy levels, pho¬ 


toionization rates, recombination rates, and Einstein co¬ 
efficients were calculated using the Flexible Atomic Code 
(FAC) 0.7.6 ( 1^ pOO^ . Some He-like emission line 

wavelengths were replaced by more accurate values from 
Verner et al. (1996b). 

We use a more approximate treatment of Fe xvii-xxiv. 
Radiative recombination (RR) and dielectronic recombi¬ 
nation (DR) line power for these ions are computed by 
(2002 ^, using atomic data from FAC. Photoionization 


and recombination are balanced to determine relative 
abundances of the i and i + 1 ionic charge states. We add 
photoexcitation from the ground level to n = 2 — 5 and 
include line emission from the resulting cascades. Similar 



Fig. 6. Composite optical Ha (red). X-ray 0.4-0.7 keV 
(green), and X-ray 0.7-1.3 keV (blue) image. The Ha im¬ 
age is from HST, while the lower resolution X-ray images 
are from Chandra. Higher ionization X-ray emission line 
regions are enveloped by a cocoon of lower ionization H 
II filaments. Regions of strong Ha emission in the NLR 
correspond to the softer X-ray emission regions in Fig. ^ 
H II regions outside the photoionization cone delineate the 
spiral arms in the starburst disk of the host galaxy. 


to the H-like and He-like ions, the relative fluxes of pho¬ 
toexcitation and recombination lines depend strongly on 
the column density and can therefore be used to constrain 
it (Sako et al. 2002). 

We start our fit by fixing the Doppler h parameters in 
the nuclear region at the values in Table ^ and fixing 
the emission redshift at its mean value < z >= 0.0025. 
The Ne x, S xv, and Fe xvii-xxiv lines are all strongly 
blended. We fix their widths at 6 ^ 800 km s“^, which 
appears to give a good overall model fit. Because the Fe 
XXV lines are strongly blended and the n > 2 Fe xxvi 
lines are not detected, the column densities of these ions 
are not constrained and we do not include them in our 
fit. We fit Gaussian profiles for the the S i and Si i Ka 
lines, so that they don’t skew the results for lines they 
are blended with. A fixed scattered continuum component 
(F = 1.8) is added to the model, extrapolated from the fits 
of the HEG continuum (see below), so that the emission 
line fluxes aren’t overestimated. 

Next we fit the RRC profiles to determine the elec¬ 
tron temperature. (The width of an RRC in eV gives a 
rough estimate of the temperature of recombining elec¬ 
trons.) The O VII and Ne ix RRCs are well measured, 
with Te = 8 — 9 X 10^ K (Table A larger range of tem¬ 
peratures is probably present for other ions observed that 
have RRCs which are too weak or blended to measure. 
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Therefore, we fix Te for each of these ions at a value which 
corresponds to an ionization parameter of peak emissiv- 
ity. These temperatures were computed using the XSTAR 
photoionization code (Table ||). Finally, we ht for column 
densities and the quantity fcL^ for each ion (Table 

The fitting procedure is similar for the off-nuclear (NE 
cloud) spectrum. Since not all emission lines are resolved, 
we assume b = 800 km s“^, the mean value for the resolved 
lines. The model spectrum is convolved with a spatial pro¬ 
file with FWHM= 1.5". Because the RRCs are so weak, 
we can not derive accurate values for the electron temper¬ 
ature in the NE cloud. Again, we assume Tg values derived 
from XSTAR (Table |). 

6.2. Ionic Column Densities and Abundances 

Fig. 0 shows our photoionization model ht for the nuclear 
region. The ht is quite good, with small residuals. Notable 
exceptions are the forbidden (f) line huxes of O vii and 
Ne IX, which appear to be underestimated by about 20%. 
We conjecture that this may be due to non-Gaussian line 
shapes. In particular, the O vii lines appear to have un¬ 
derlying broad components which we do not model. The 
opposite problem occurs with O viii Ly/3, which is overes¬ 
timated because it has a narrower width than O viii Lya. 
This is difficult to understand, unless O viii Lya is broad¬ 
ened more by scattering or there are additional emission 
sources which add to the line. There appears to be un¬ 
modeled excess emission in the 8.8-9.5 A region, in the 
vicinity of the Mg xi triplet, which may affect our column 
density estimate of this ion. 

We find nuclear ionic column densities (Table in 
the range 2 x 10^® — 6 x 10^® cm“^, which are comparable 
to warm absorber column densities observed in Seyfert 1 


galaxies (e.g., Kaastra et al. 


is larger perpendicular to the axis of the ionization cone 
than parallel to it. However, it is possible that bulk mo¬ 
tions of clouds may broaden the lines without affecting 
the resonance line optical depths in individual clouds. In 
that case, our measured line widths give an upper limit to 
6, and our column densities are upper limits. As pointed 


(200C)). We do not give for- 
because it is 
40 pa- 


mal uncertainties for the column densities 
hard to calculate correlated uncertainties in a 
rameter fit. However, we do give (90%) single-parameter 
uncertainties for fcL^. We find < fcL^ >= 1.5±0.2 x 10^^ 
erg s“^ for the nuclear region. 

The large r// (> 1) values observed for Mg xi and Si 
XIII are nicely explained by the photoionization plus pho¬ 
toexcitation model. The abundances and column densities 
of these ions (Table jj) are relatively low so photoexcita¬ 
tion dominates over recombination. This is in contrast to 
the O VII and Ne ix ions, which have larger column densi¬ 
ties and saturated 2p-ls resonance lines. There is therefore 
no need to invoke collisional excitation to explain the line 
ratios. We also reemphasize the fact that the strong 3p-ls, 
4p-ls, and 5p-ls lines require photoexcitation, and would 
not be present in a collisionally excited plasma. 

The column densities we derive for the nuclear re¬ 
gion are a factor of a few higher than those derived by 
Kinkhabwala et al. ( 2002| ), who use a smaller radial tur¬ 
bulent velocity of = 200 km s“^ {b = 280 km s“^) and 
an additional line of sight velocity broadening of = 400 
km s“^. It seems unlikely to us that the turbulent velocity 


out by 3ako et al. (2002), since the line fluxes are pro¬ 
portional to column density, fit column densities will be 
biased toward large values. 

Fig. H shows our photoionization model ht for the NE 
cloud. The model does a fairly good job of representing 
the data, except for residuals in the 11.7-12.4 and 15-15.2 
A regions. The Fe xxi emission falls in the former band, 
precluding an accurate determination of the column den¬ 
sity. The Fe xvii column density is underestimated by 
the model {N[ = 4 x 10^® cm“^). This is apparent from 
comparing the strengths of the 3s-2p (17.1 A ) and 3d-2p 
(15.0 A ) lines which indicate a more consistent value of 
Ni = 2 X 10^^ cm“^. IMP tried to decrease the unidenti- 
Hed residual at 15 A by strengthening the 3d-2p line with 
respect to the 3s-2p line, leading to a column density es¬ 
timate which is too low. The Al Xii, S xv, and Fe xxiv 
emission are too weak to constrain column densities. We 
also hnd an unexplained shift in the wavelength of the 
O VII f line, though the other lines from this series are 
consistent with zero redshift. 

The covering fraction of the NE cloud is indistinguish¬ 
able from the nuclear region, with /cLx = 9.0 ±0.7 x 10^® 
erg s“^. The small dispersion (30-50%) in fcL^ in both 
regions indicates that all ions have a similar covering frac¬ 
tion, assuming the source is isotropic within the ionization 
cone. 

We estimate the equivalent hydrogen column iVn and 
emission measure EM= J nldV for each ion (Table 
assuming solar abundances and ion fractions modeled by 
XSTAR. The mean column densities averaged over all ions 
are A^h = 5 ± 1 x 10^^ cm“^ in the nuclear region and 
iVn = 1.4 ± 0.3 x 10^^ cm“^ in the NE cloud. The factor 
of 4 lower mean column density in the NE cloud accounts 
for the weaker forbidden lines and RRCs relative to the 
nuclear spectrum. Photoexcitation is relatively stronger 
because the resonance lines are not saturated. The relative 
column densities of He-like Mg and Si are even lower in 
the NE cloud, accounting for their very large r/f ratios. 
Such large ratios are impossible to produce by collisional 
ionization in a low-density plasma. In contrast, relatively 
large Fe ion column densities in the NE cloud account for 
the strong Fe L-shell emission. 

Relative elemental abundances (Table ||) are esti¬ 
mated by summing the ionic columns and taking the ra¬ 
tio with respect to neon. These are rough estimates since 
we are comparing clouds with a large range of ionization 
states, and ignore possible abundance gradients. Except 
for oxygen and iron, the relative abundance patterns for 
both the nuclear region and NE cloud are similar (within 
a factor of 2) to solar. Oxygen appears to be under- 
abundant by a factor of 4.0 in the nuclear region, but is 
near solar abundance in the NE cloud. Oxygen depletion 
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Table 4. NGC 1068 relative covering fraction, column 
density, and emission measure of the nuclear region (Nuc.) 
and NE cloud (NE). The quantity fcL^ is normalized to 
the mean for the nuclear region, < fcLx >= 1-5 x 10^^. It 
therefore gives the relative covering fraction, assuming Lx 
is the same for all ions. Ionic column densities Ni are given 
in units of 10^^ cm“^ and equivalent hydrogen column 
densities in units of 10^^ cm“^. As explained in the text, 
IMP underestimates the column density of Fe XVII in the 
NE cloud. The value presented in parentheses here gives a 
better match to the spectrum. The emission measure EM 
is in units of 10®^ cm“^. 


Ion 

Nuc. /cLx 

Ni 

Nh 

EM 

NE /cLx 

Ni 

Nu 

EM 

N VII 

0.8±0.2 

49 

8.3 

16 

0.7±0.2 

27 

4.6 

11 

O VII 

1.7±0.2 

37 

0.77 

20 

0.9±0.2 

30 

0.62 

7.6 

O VIII 

0.51±0.09 

59 

1.3 

2.1 

o 

d 

d 

o 

d 

98 

2.2 

1.8 

Ne IX 

0.45±0.07 

49 

7.4 

7.2 

0.3±0.1 

14 

2.1 

1.4 

Ne X 

1.1±0.2 

6.1 

1.0 

0.76 

0.6±0.2 

5.3 

0.87 

0.36 

Mg XI 

0.39±0.07 

23 

12 

3.9 

0.8±0.3 

1.2 

0.61 

0.62 

Mg XII 

0.62±0.09 

10 

5.4 

1.4 

0.9±0.2 

1.7 

0.91 

0.33 

Al XII 

1.3±1.0 

0.20 

1.4 

1.4 





Si XIII 

0.8±0.1 

19 

11 

3.6 

0.8±0.3 

2.3 

1.3 

0.42 

Si XIV 

0.7±0.1 

14 

8.3 

1.5 

0.6±0.3 

2.7 

1.6 

0.22 

S XV 

0.6±0.4 

20 

25 

3.4 





Fe XVII 

2.5±0.6 

0.44 

0.30 

0.72 


(2) 



Fe XVIII 

1.5±0.4 

0.81 

0.66 

1.0 

0.4±0.1 

1.5 

1.2 

0.51 

Fe XIX 

0.9±0.3 

1.1 

2.6 

0.56 

0.8±0.2 

0.98 

2.3 

0.43 

Fe XX 

0.8±0.3 

1.3 

1.0 

0.27 

0.5±0.2 

0.83 

0.66 

0.11 

Fe XXI 

0.7±0.3 

2.3 

3.2 

0.15 





Fe XXII 

1.9±0.7 

0.49 

0.35 

0.08 

0.9±0.6 

0.38 

0.28 

0.027 

Fe XXIII 

0.8±0.5 

1.1 

0.97 

0.09 

0.4±0.4 

0.81 

0.71 

0.031 

Fe XXIV 

1.0±0.9 

0.62 

0.89 

0.05 






Table 5. NGC 1068 relative elemental abundances. 


Element 

Nuclear 

Region 

NE cloud 

Solar 

N 

> 0.89 

> 1.40 

0.91 

0 

1.74 

6.63 

6.92 

Ne 

1.00 

1.00 

1.00 

Mg 

0.60 

0.15 

0.31 

Al 

> 0.004 


0.024 

Si 

0.60 

0.26 

0.29 

S 

> 0.36 


0.13 

Ee 

0.15 

> 0.34 

0.38 


by a factor of 5 was previously inferred (Marshall et al. 
1993|) based on the weakness of the O viii Lya line. A 


similar depletion by a factor of 4.0 has been inferred for 
nuclear molecular clouds based on a large HCN/CO in¬ 


tensity ratio (Sternberg et al. 1994). The Fe abundance 


is depleted by a factor of 2.5 in the nuclear region, but is 
near solar abundance in the NE cloud. 


6.3. Limits on Collisionally Ionized Plasma 


Our observations confirm that line emission from the X- 
ray NLR in NGC 1068 is dominated by warm (Te = 
8 X 10"*) photoionized gas. However, past observations in¬ 
dicate that shock-heating may also be important in the 
NLR. The NE radio lobe terminates at a distance of ~ 5" 
from the nucleus and has a sh arp leading edge, possibl y 
corresponding to a shock front ( Wilson fc Ulvestad 1987 ). 
There are also indic ations that the NLR clouds decelerate 
on the 2 — 4" scale ( Crenshaw fc Kraemer ^000 ), perhaps 
by plowing through the ISM of the host galaxy. Other 
studies have shown that photoionization by emi ssion from 
shocks may be im portant in the optical NLR ( Dopita fc 
Sutherland 1995). In addition, pressure confinement of 
the NLR clouds may require the presence of a hot com¬ 
ponent with Te cs 10^ K ( Krolik et al. 1981 ; Ogle et al. 

Spool) . 


Adding a collisionally ionized (Cl) component to the 
spectrum would increase the strength of the 2p-ls (a) lines 
of H-like and He-like ions. The contribution to forbidden, 
inter-combination, and high-order Lyman lines would be 


relatively weaker (Kinkhabwala et al. 2002|). Since the 


line ratios observed in NGC 1068 are closely consistent 
with a photoionized plasma, the 2tT uncertainties in the 
n = 2 — 1 line fluxes give reasonable upper limits to the 
contribution from CL (Similar estimates could be made 
from the residual broad O vii or O viii emission.) We 


use the emissivities and peak temperatures of Mewe et al 
(|1985| ), who assume solar abundances. From the O vii 
2p-ls line (Table H), the 2cr upper limit for the emission 
measure of a collisionally ionized component with T = 
2 X 10® K in the nuclear region is EM< 6 x 10®^ cm“®. 
This is 30 times smaller than the emission measure of the 
photoionized component (Table ^). The upper limit for a 
hotter collisionally ionized component with T = 1 x 10^ K 
is less restrictive and comparable to the emission measure 
of the photoionized Mg xii Lya line: EM < 3 x 10®® cm“®. 

The corresponding limits on Cl plasma for the NE 
cloud are EM< 5 x 10®^ and EM< 1 x 10®® for T = 
2 X 10®, 1 X 10^, respectively. These are smaller by factors 
of 15 and 2.5, respectively, than the observed emission 
measures of photoionized O Vii 2p-ls and Mg Xii Lya. 
A Cl component would also emit thermal bremsstrahlung 
continuum and strong Fe L lines ( ^ako et aT7 2000). These 
are more difficult to constrain. 

While we rule out collisional ionization as a major 
contributor to the spectrum, a weak component could 
be present. As we noted above, there is a region in be¬ 
tween the nuclear region and NE cloud (Fig. ^ with a 
hard spectral color, similar to that in the starburst region. 
This region is relatively free of NLR clouds in the HST [O 
III] image (Macchetto et al. 1994). With a higher S/N 


Chandra grating observation at the appropriate roll an¬ 
gle, it should be possible to get a spectrum of this region 
to see if it contains any hot, collisionally ionized plasma. 
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Fig. 7. Photoionization model of MEG spectrum of nuclear region. Ion spectra are fit with IMP, to determine column 
densities and temperatures (see text). The ratios of photoexcitation to recombination lines are the primary diagnostics 
for column density. The widths of the RRCs are a measure of the electron temperature. 
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Fig. 8. Photoionization model of MEG spectrum of the NE cloud. Similar to Fig. 
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6.4. Characterization of the NLR Outflow 


The results of the last two sections show that the X-ray 
NLR clouds are ionized and heated primarily by X-ray 
photons. The luminosity of the photoionizing source is 
consistent with the hidden AGN (Sect. 6.6). There is also 
no indication that the radio jet or lobe is a significant 
source of X-rays, otherwise they would appear directly 
in our X-ray images and spectra. Mechanical heating and 
ionization by shocks are apparently weak, if at all present. 
Any contribution to the nuclear X-ray spectrum from col- 
lisionally ionized gas, either from a starburst or shocked 
by the radio jet has yet to be detected. 

On the other hand, there is abundant morphological 
and kinematic evidence that the radio jet plays a role in 
shaping the NLR and that it influences cloud dynamics. 
This leads us to believe that the radio jet interaction with 
the NLR is a rather gentle process and indicates a subsonic 
jet. An analogy may be drawn to the radio jets which 
evacu ate bubbles in the I GM surrounding radio galaxies 
(e.g., Fabian et al. ( 2000 )). 


Our results also indicate that radiation power dom¬ 
inates over kinetic power output in NGG 1068. Kinetic 
power from the radio jet can at most provide a small frac¬ 
tion of the heating. However, radiation pressure and jet 
pressure may still be comparable in affecting the dynam¬ 
ics of the NLR. The kinem atics of the NLR suggest two 
different regions of influence. Crenshaw fc Kraemer ( ^000 ) 
model the inner 1.5" of the optical [O ill] emission line re¬ 
gion with radially accelerating clouds in a hollow bi-cone. 
This structure may have naturally formed by ablation and 
radiative acceleration of matter from the torus. In the re¬ 
gion of the NE cloud, the radio lobe may be expanding 


into the host galaxy disk (Gecil et al. 2001). 

We use a grid of XSTAR 2.1 (Kallman 2001) pho¬ 
toionization models covering the parameter space log ^ = 
0.0 — 4.0 (^ = L/riei'^) to estimate temperatures, ion frac¬ 
tions, and line emissivities in the X-ray NLR. We assume 
sol ar abundances. The in put SED is similar to that used 
by Krolik fc Kriss ( 2001 ). The X-ray spectrum is a power 
law with photon index T = 1.8 and high energy cut-off at 
100 keV. We assume that each emission line is produced 
in a region of peak emissivity to derive an ionization pa¬ 
rameter and temperature for each ion. 

Fig. ||a shows the distribution of A^h in the nuclear 
region (from Table jj) as a function of f and T. There is a 
large amount of scatter, which we suggest is primarily due 
to non-solar abundances. In particular, relatively low oxy¬ 
gen and iron column densities suggest that these elements 
are under-abundant. The distribution covers a large range 
in ionization parameter, corresponding to a factor of 80 in 
density or a factor of 9 in distance from the central source. 

We compute the differential emission measure (DEM) 
distribution (Fig. ^) from the EM values (Table ^ 
assuming d{logf) = 1.0,0.5,0.25 for H-like, He-like, and 
Fe L ions, respectively. These d{logf) are estimated from 
the widths of the i 1 ion fraction curves fi+iif) com¬ 
puted with XSTAR. Assuming equal-size errors for all 


the data points, we find a best-fit slope of -1.0 and an 
intercept of log{d{EM)/d{logf)) = 65.6. The slope can 
be understood if we approximate the A^h distribution by 
.^h(C) ~ constant. Then the emission measure drops with 
ionization parameter, following EM^ (Fig- ^)- Using 
the dehnitions of emission measure and ionization param¬ 
eter, this also implies that the product of the mass and 
X-ray flux is constant, 

-^(C)F'x rieVFx ^ f x EM = constant. (1) 

In other words, at any given distance from the source, 
there are equal masses of plasma at all observed ionization 
parameters. This implies a large range of density at any 
given location 
( ^OOlD suggest 

ture may result from marginal stability along the vertical 
branch of the cooling curve. It is also apparent from the 
soft X-ray image (Fig. ^) that the emission from the X- 
ray NLR is clumpy and azimuthally asymmetric, with a 
range of colors which isn’t a simple function of distance 
from the nucleus. 


( Kinkhabwala et al. 2002|) . Krolik fc Kriss 


that a large range in density and tempera- 


6 . 5 . The Torus 


The obscuring region in NGG 1068 can be pictured as a 
collection of dusty molecular clouds surrounding the nu¬ 
cleus. A conical illumination pattern for the ionization 
cone is most simply created by a toroidal ring, but more 
complicated geometries may produce a similar effect. The 
200 X 160 pc (projected) central region of the Ghandra 
image (Fig. ^ shows 3-8 keV emission from clouds sur¬ 
rounding the hard X-ray peak. This emission may be re¬ 
flected continuum and Fe Ka fluorescence. The lack of 
1.3-3 keV line emission from this region is consistent with 
a low ionization parameter and high density. 

The extended hard X-ray emission surrounding the nu- 
cleus coincides wi t h the molecular CO ring measured by 
^chinnerer et al.' ( 2000| ). The two brightest CO spots in 
this ring are located east and southwest of t he nucleus 
and contain a n estimated mass of 5 x 10® Mq. ^chinnera 
let al. (|200C| ) model the kinematics of the molecular gas 
with a highly warped disk. The axis of the disk is 70° 
to the line of sight and has a PA of 20° at a distance 
of 1" from the nucleus. A warped disk geometry can cre¬ 
ate a highly irregular X-ray illumination pattern, perhaps 
leading to the observed irregularity in the hard X-ray ring. 
Alternatively, the blobs in the ring may correspond to dis¬ 
crete giant molecular clouds. Much closer to the nucleus, 
H 2 O masers outline an edge-on portion of the molecular 
disk with an outer radius of 30 mas (2pc) and PA= 4°, 
and perhaps the inner surface of the torus ( pallimore et 
|al. 1996; Greenhill et al. 1996). It is likely that the col¬ 
umn density of the inner disk is Compton-thick, sufficient 
to completely obscure the direct line of sight to the nu¬ 
cleus. However, the molecular clouds at 100 pc scale also 
appear to be important in defining the geometry of the 
ionization cone. 
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with a PEXRAV ( Magdziarz fc Zdziarski 1995|) com- 


cm 


keV ^ and 4±lxl0 "^phs ^ cm ^ keV 


nosity from the ionized region is r/cLx = 4 ± 1 x 
erg s“^. 

The equivalent width of the Fe i Ka line is 1.5 ± 0.3 
keV, consistent with the value measured by BBXRT 


( Marshall et al. 1993 ). Such a large equivalent width is 
in agreement with Monte Carlo simulations of neutral re¬ 
flection from the inner wall of a Compton-thick molecular 
torus ( IChisellini et al. 1994 ). It is inconsistent with the 
smaller equivalent width of ~ 0.5 keV p redicted for reflec- 
tion from the ionized scattering region (Krolik & Kallman 


1987). 


6.6. Ionized Scattering Region 

An ionized scattering region which gives an indirect view 
of the hidden nucleus is a necessary component of the 
Seyfert unification theory. Ultraviolet light scattered from 
the inner 4" of NGC 1068 has a wavelength-independent 


polarization (Antonucci et al. 1994). This is strong evi¬ 


dence for electron scattering as opposed to dust scattering. 
The optical scattering fraction is estimated at r/c = 0.015 


from spectropolarimetry of the H/3 line (Miller et al. 


Fig. 9. a) Column density and b) differential emis¬ 
sion measure (DEM) distributions in the nuclear region. 
Ionization parameter ^ and temperature are from XSTAR 
photoionization models, assuming maximum emissivity at 
solar abundance. The distributions cover a large range (2 
decades) of ionization, with column densities of similar 
magnitude in each ionization state. The relatively low col¬ 
umn densities of oxygen (X’s) and iron (squares) suggest 
they are under-abundant. The best fit model to the DEM 
distribution gives EM ~ indicating equal mass at 
each ionization parameter. 


A further confirmation of the obscuration model is the 
reflection spectrum from the nuclear region (Fig. We 
first fit the 1.2-5 A (2.5-10 keV) hard X-ra y cont inuum 


1991). The electron-scattering regio n is thought to be asso - 
ciated with the ionized X-ray NLR (|Krolik fc Kriss 1995 ). 


ponent representing neutral reflection from the molecular 
torus. We also include Gaussian emission lines at the en¬ 
ergies of Fe I (Kof and K/3) and Fe XXV (r and f), and 
account for Galactic absorption. We use the C-statistic 
since there are few counts per bin in the HEG spectrum. 
We fit a neutral reflection component with an incident 
power law index of T = 2.8 ± 0.2. This unusually steep in¬ 
dex and significant residuals in the 4-5 A region indicate 
that an extra continuum component is present. Next we 
add a power law component from ionized reflection. The 
incident power law index is fixed to T = 1.8 for both re¬ 
flection components, since it is not well constrained by the 
data. The best fit normalizations at 1 keV for the neutral 
and ionized reflection components are 1.2 ± 0.2 x 10 “^ ph 


As shown in the last section, there is evidence for electron 
scattering from this region in the hard X-ray band. The 
scattered X-ray continuum and emission line flux from the 
nuclear region give estimates of t/cLx and /cTx- From the 
ratio of these, we find an electron scattering optical depth 
of r = 0.27 ± 0.08 (90% confidence.) Gombining this re¬ 
sult with the optical scattering fraction, the source cov¬ 
ering fraction in the nuclear region is fc — 0.06, and the 
luminosity of the hidden nucleus is ~ 2.7 x 10"*^ erg 
s“^. The mean equivalent Thomson scattering depth as¬ 
sociated with any given ion , r = 3.2 x 10“^ (derived from 
< Vh >, Table ^), can account for 12% of the scattered 
X-ray continuum. Integrating over all ionization states, we 
can account for the total electron scattering column (Fig. 

©• 

The polarized broad H/3 line from the central 4" of 
NGC 1068 appears to be broadened by the thermal motion 
of scattering electrons with mean temperature < >= 

6.7 x 10^ K (Miller et al. 1991). Note that the higher 
value < Te >= 3 x 10® K which is often quoted assumes 
a monochromatic incident line, whereas the lower value 
we quote here is for a broad incident line width estimated 
from the dust scattered line in the off-nuclear knot. The 
electron temperature we measure for the O vii ion (Tg = 
8^2 X lO"* K) is in good agreement with the value from 
optical spectropolarimetry. The ^ 400 km s“^ redshift 
of the polarized broad H/3 line, attributed to the outflow 


of the electron scattering region (Miller et al. 1991), is 


respectively. The corresponding scattered 2-10 keV lumi- 


consistent with the mean blue shift of the X-ray emission 
lines, < V >= —410 km s“^. 

A multi-temperature electron scattering region with 
logT(Ar) = 4.5 — 6.0 is indicated by the observed column 
density distribution (Fig. We integrate over model 
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temperature distributions (Fig. 0a) to predict the op¬ 
tical H/3 profile scattered from the X-ray NLR. A simple 
pedestal function approximates a flat temperature distri¬ 
bution with a high temperature cut-off Tq and low tem¬ 
perature cutoff at logT(Ar) = 4.5. The area under the 
curve is normalized to give r = 0.27. We try two different 
models, with logTc = 5.0, 5.95. We integrate the ionic col¬ 
umn densities over the model temperature distributions, 
and divide by the ionic abundances from our XSTAR mod¬ 
els to find column density distributions (Fig. 0b) which 
are consistent with the observed distribution (Fig. 

We use a Gaussian incident BLR profile with FWHM= 


A-- - 


-^-^ 


2900 km s ^ (Miller et al. 1991) to compute the scat¬ 


tered H/3 profile. We also include turbulent broadening of 
h = 1000 km s“^, the mean value for the X-ray emission 
lines (Table ^). The resulting profiles are shown in Fig. 
0c, and compared to the observed polarized broad H/3 


line width FWHM= 4400 km s ^ (Miller et al. 1991). 



The model with logTc = 5.0 gives a better match to the 
observed H/3 profile than logTc = 5.95. This indicates 
that there is relatively little plasma at hotter tempera¬ 
tures, or it would broaden the H/3 line more than is ob¬ 
served. Broadening by turbulent motions in the NLR has 
an effect on the profile width for logTc = 5.0, but is in¬ 
significant for logTc = 5.95 (Fig. [l^,d). 

One concern is that the optical-UV bremsstrahlung 
from the ionized scattering region not exceed the observed 
continuum in that wave band. Miller et al. (1991) found 



that this constraint is satisfied in their models, which as¬ 
sume a smooth distribution of matter. However, if the 
scattering clouds are clumpy, then bremsstrahlung is en- 


hanced by a factor of < > / < Ue (|Krolik fc 


Kriss ||;995[)r We can make a direct estimate of the optical- 

UV bremsstrahlung by integrating over the DEM(T) curve 
(Fig. ib). We predict peak emission of vF^ = 8.4 x 10 
erg s~^ cm~^ at A = 300 0 A. This is 40 times lower than 
the [Miller et al. ( 1991 ) estimate of < 3.4 x 10 


Fig. 10. a) Model temperature distributions of optical 
depth in the nuclear scattering region, assuming two differ¬ 
ent cut-off temperatures (logTc(AT) =5.0,5.95; dash and 
dot-dash), b) Corresponding column density distributions, 
both consistent with the observed distribution in Fig. 9. c) 
Scattered H/3 profiles, integrated over the two model tem¬ 
perature distributions in panel a. The observed (incident 
and scattered) line profiles are the solid curves (Miller et 


al. 1991). The dashed and and dot-dashed curves show 


erg s“^ cm“^ (970A), which assumed a radial density dis¬ 
tribution of rZe ~ r~^. Hence bremsstrahlung emission is 
much weaker than the scattered optical-UV light and will 
be very difficult to detect. 

Our Chandra observations confirm the hypothesis that 
the X-ray NLR and the optical/X-ray electron scattering 
region are one and the same. Both occupy the same ex¬ 
tended region in the central 3" of the galactic nucleus, and 
have roughly the same optical depth and electron tem¬ 
perature. The expected variability timescale is the light¬ 
crossing time, which is 240 years. This is apparently at 
odds with a recent report of 20% variability in 4 months 
of the scattered X-ray continuum and 60% variability of 
the Fe xxv (6.7 keV) emission line flux (Colbert et al. 
2002| ). However, we can not rule out the possibility that 


the continuum variability comes from a point source which 
contributes a significant fraction of the flux. There are 
not enough photons in our Chandra observation to obtain 
separate images of the Compton-reflected and electron- 
scattered components of the X-ray continuum. Similarly, 
we can not spatially distinguish the regions of Fe i and Fe 


the predicted scattered profiles, including both thermal 
and turbulent broadening (1000 km s“^). d) Similar to c), 
without turbulent broadening. 


xxv Ka emission. It will be important to follow up with 
additional high-resolution Chandra spectroscopy to see if 
the Fe xxv line variability is confirmed. 

7. Conclusions 

Chandra HETGS X-ray spectroscopy and imaging of NGC 
1068 provide new diagnostics of the spatially resolved 
narrow-line region and support the Seyfert unification 
model. The X-ray spectrum contains emission lines from 
a large range of ions in the NLR, and fluorescence from 
the molecular torus. The X-ray emitting plasma in the 
extended NLR is primarily photoionized and there is no 
evidence for a collisionally ionized component. Our previ¬ 
ous claim of such a component in NGC 4151 (jOgle et al. 


^000 ) was confused by the effects of photoexcitation. 

Observations of the narrow Fe Ka line show that it 
is produced predominately in the nuclear region, near the 
location of the hidden nucleus. Its equivalent width is con¬ 
sistent with obscuring torus models. It also appears that 
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we have detected and resolved the outer regions of the 
molecular torus in hard X-rays. The hard X-ray continuum 
spectrum from the nuclear region consists of a Compton 
hump from neutral scattering and an electron scattered 


compo rent. 


We derive ionic column densities, equivalent hydrogen 
column densities, and emission measure distributions for 
both the nuclear region and a cloud 3" NE of the nucleus. 
The column density of the X-ray NLR can account for the 
entire electron scattering column in the nuclear region. 
In addition, its temperature and velocity are sufficient to 
cause the observed broadening and redshift of the scat¬ 
tered optical broad lines. These facts are all consistent 
with electron scattering in the X-ray NLR providing the 
view of the hidden Seyfert 1 nucleus. 

The results for the X-ray NLR in NGC 1068 may also 
have important implications for the properties of ionized 
(warm) absorbers seen in the spectra of many Seyfert 1 
galaxies. Similar velocities and ionization parameters are 
found in warm absorbers and the X-ray NLR. We see 
from Chandra observations of NGC 1068 that AGN out¬ 
flows can be clumpy and show a wide range of ionization 
log 5 = 1.0 — 3.0. The differential emission measure dis¬ 
tribution indicates roughly equal mass in each ionization 
state in this range. The emission seen in various UV and 
X-ray lines, though occupying roughly the same region, 
must come from clouds of greatly different ionization, tem¬ 
perature, and density. It is important to understand the 
equilibrium among the different phases to develop a com¬ 
plete picture of AGN outflows. 
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